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The dynamical spin susceptibility is studied in the magnetically-disordered phase of heavy- 
Fermion systems near the antiferromagnetic quantum phase transition. In the framework of the S = 
1/2 Kondo lattice model, we introduce a perturbative expansion treating the spin and Kondo-like 
degrees of freedom on an equal footing. The dynamical spin susceptibility displays a two-component 
behaviour in agreement with the inelastic neutron scattering (INS) experiments performed in CeCu§, 
Cei- x La x Rii2Si2 and UPt^: a quasielastic q-independent peak as in a Fermi liquid theory, and a 
strongly q-dependent inelastic peak typical of a non-Fermi liquid behaviour. Very strikingly, the 
position of the inelastic peak is found to be pushed to zero at the antiferromagnetic transition. The 
picture is consistent with the neutron cross sections observed in INS experiments. 

PACS numbers: 71.27.+a, 75.20.Hr, 75.40.Gb 



I. INTRODUCTION 



One of the most striking properties of heavy Fermion compounds discovered these last years is the existence 
of a quantum phase transition driven by composition change (at xc = 0.1 in CeCu§- x Au x and xc = 0.08 
in Cei- x La x Ru2Si2), pressure or magnetic field. It has been largely discussed in various theoretical approaches 
II^J . Important insight is provided by the evolution of the low temperature neutron cross section measured by 
Inelastic Neutron Scattering (INS) experiments when getting closer to the magnetic instability. The experiments 
performed in pure compounds CeCue and CeRu2S%2 by Regnault et al |^] and Aeppli et al p| have shown the 
presence of two distinct contributions to the dynamic magnetic structure factor: a q-independent quasielastic 
component, and a strongly q-dependent inelastic peak with a maximum at the value w max of the frequency. 
The former corresponds to localized excitations of Kondo-type while the latter peaked at some wavevector Q is 
believed to be associated with intersite magnetic correlations due to RKKY interactions. The frequency-width of 
the quasielastic and inelastic peaks respectively define the single-site and intersite relaxation rate T ss & n d T/s- 
Such features have also been observed in UPt$ and called as "slow" and "fast" components by Bcrnhoeft and 
Lonzarich 0. Later on, INS experiments have been performed with varying compositions as in C e\- x La x Ru2S%2 
M . It has been observed a narrowing of the single-site relaxation rate Tss when getting closer to the magnetic 
transition. At the same time, both the position of the inelatic peak w max and the intersite relaxation rate Tjg 
drastically decrease when getting near the magnetic instability. Table 1 reports the values of TsSj Frs and w min 
for the different compounds. 

Any theory aimed to describe the quantum critical phenomena in heavy-Fermion compounds should account 
for the so-quoted behaviour of the dynamical spin susceptibility. We start from the Kondo lattice model which 
is believed to describe the physics of these systems. We refer to the recent paper of Tsunetsugu et al M for a 
review of the model. As already pointed out by Doniach in his initial paper Jlfl], the main features result of 
the competition between the Kondo effect and the RKKY interactions among spins mediated by the conduction 
electrons. Most of the theories developed so far JTT|-p^| agree with the existence of a hybridization gap which splits 
the Abrikosov-Suhl or Kondo resonance formed at the Fermi level. The role of the interband transitions has been 
outlined for long in order to explain the inelatic component of the dynamical spin susceptibility. For instance, the 
theories based on a 1 /N expansion |l6|-|2C|] (where N is simultaneously the degeneracy of the conduction electrons 
and of the spin channels) predict a maximum of x {kp,ui)/u> at u max of the order of the indirect hybridization 
gap pi] ]. However, the 1/N expansion theories present serious drawbacks: (i) the spin fluctuation effects are 
automatically ruled out since the RKKY interactions only occur at the following order in 1 /N 2 E^j, (ii) they then 
fail to describe any magnetic instability and hence the quantum critical phenomena mentioned above and (iii) the 
predictions for a-> max and the associated relaxation rate cannot account for the experimental observations near the 
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magnetic instability. An improvement brought by Doniach |2j| consists to consider the 1/N 2 corrections in an 
instantaneous approximation: it gives back the ladder diagram contribution to the dynamical spin susceptibility 
and then accounts for the spin fluctuation effects. Other approaches were proposed in Ref. [^4] ^(| . But still the 
predictions for the frequency dependence of the dynamic magnetic structure factor presents a gap of the order 
of the hybridization gap whatever the value of the interaction is. On the other hand, in front of the difficulties 
encountered when starting from microscopic descriptions, various phenomenological models (as the duality model 
of Kuramoto and Miyake [^7| and reference JtJ) have been introduced to describe both the spin fluctuation and 
the itinerant electron aspects with some successful predictions as the weak antiferromagnetism of these systems. 

In this paper, we develop a systematic approach to the Kondo lattice model for S = 1/2 (N — 2) in which the 
Kondo-likc and the spin degrees of freedom are treated on an equal footing. The presented approach shows some 
similarities with earlier works P|,ff5| . But while Ref. 0,^5| essentially describe the phase diagram of the Kondo 
lattice at a mean-field level, we focus on the effects of spin fluctuations in the magnetically-disordered phase hence 
bringing the spin- fluctuation and the Kondo effect theories together. The saddle-point results and the gaussian 
fluctuations in the charge channel are consistent with the standard 1/N theories. In addition, the gaussian 
fluctuations in the spin channel restore the spin fluctuation effects which were missing in the 1/N expansion. 
The general expression of the dynamical spin susceptibility that we derive reproduces some of the features 
postulated in the phenomenological models. It presents a two-component behaviour: a quasielatic component 
superimposed on an inelatic peak with renormalized values of the relaxation rates, susceptibilities and w max - In 
a very striking way, w max is pushed to zero and the inelastic mode becomes soft at the antiferromagnetic phase 
transition with vanishing relaxation rate. Predictions are quantitatively compared with experimental results. 
The quasielastic peak is typical of a Fermi liquid while the other mode breaks the Fermi liquid description. 
Our approach might offer new prospects for the study of the quantum critical phenomena in the vicinity of the 
antiferromagnetic phase transition. 



II. PRESENTATION OF THE APPROACH 



We consider the Kondo lattice model (KLM) constituted by a periodic array of Kondo impurities with an 
average number of conduction electrons per site n c < 1. In the grand canonical ensemble, the hamiltonian is 
written as 

H = £kc{ a c ka + Si ■ Y Aa T °° ,c i<r' - M^s( c L c k<? ~ n c ) (1) 

ka i aa' ka 

in which Si represents the spin (S = 1/2) of the impurities distributed on the sites (in number Ns) of a 
periodic lattice; ct CT is the creation operator of the conduction electron of momentum k, spin quantum number a 
characterized by the energy tk — — X) ex P (*k.Ry ) and the chemical potential fi; t are the Pauli matrices 

<i,j> 

{t x , t v ', t z ) and r° the unit matrix; J is the antiferromagnetic Kondo interaction (J > 0). 
We use the Abrikosov pseudo-fermion representation of the spin S^: S,; = ^ f\ a T aa 1 fia 1 ■ The projection into 

aa' 

the physical subspace is implemented by a local constraint 

ia 

A Lagrange multiplier Xi is introduced to enforce the local constraint Qi. Since [Qi,H] = 0, Xi is time- 
independent. 

In this representation, the partition function of the KLM can be expressed as a functional integral over the 
coherent states of the fcrmion fields 



VdaVfiadXi exp 



C{r)dT 



(3) 



where the Lagrangian C(t) is given by 
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C(t) = C (r) + H (t) + Hj(t) 



Hq(t) = Y e k c l a c kl7 - fiN s ( — Y c\ a c k(J - n c j + Y X iQi 

kcr \ ka / i 



Hj{t) = jJ2 S fi' S c 



with S Ci = 4o T <j<J lC i<y' and S/i = Si 

We perform a Hubbard-Stratonovich transformation on the Kondo interaction term Hj[t). Since more than 
one field is implied in the transformation, an uncertainty is left on the way of decoupling. We propose to remove 
it in the following way. First, we note that Hj(t) may also be written as 



3j ^ j 

Hj{t) = — ~Y2 n fc,n c f, + -z Y S / C i ' Sc /> 



(4) 



where S/ Ci = fi a T <?<?' c i<?' anc ^ n fci — fla T aa' c i<y' (respectively S c f t and n c / ; their hermitian conjugate). 



The Kondo interaction term is then given by any linear combination of J^ S/i • S C i (with a weighting factor 

i 

x) and of the term appearing in the right-hand side of Equation (^) (with a weighting factor (1-x)). x is chosen 
so as to recover the usual results obtained within the slave-boson theories. One can check that this is the case 
for x = 1/3. The Kondo interaction term is then given by 



Hj(t) = JsY ( S /< ' S c« + S /c, ' S c/J ~ Jc Y n fc* n cf> 



(5) 



with Js = J/4 and Jc = J/3. 

Performing a generalized Hubbard-Stratonovich transformation on the partition function Z makes the fields 
$i, (for charge) and £ c . appear (omitting the fields associated to S/ Ci , S c /J. We get 



with 



Z = J d$ i d$*d$ f .d£ c .Vc i<T Vf i(T d\ i exv - J C'{r)dr 



£{t)=C {t)+H q {t) + H' j {t) 



(6) 



ajW " ^ [ Ci ° h ° ' \ Jc^tL, -Js** ■ t 



fio 



JcY^i+JsY^-^ 



A. Saddle-Point 



The saddle-point solution is obtained for space and time independent fields $o > Ao, £/ ancl £c - I n the 
magnetically-disordered regime (£^ = £ co = 0), it leads to renormalized bands a and (3 as schematized in Figure 
1. Noting ctq*^ = Jc^o*^ and £/ = Ao, ot\ a \§ > and /3| CT |0 > are the eigenstates of 
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_ ( 3 T + £k CTq 



G -(k,r)=p— (7) 

with respectively the eigenenergies (d T + E^) and (d T + E£). In the notations: x k = e fe — ey, = — Ef 
and Afe = y/ 'x\ + 4er 2 , we get 

y± - (x k ± A fe ) /2 (8) 

Let us note [/^ the matrix transforming the initial basis {c\ a fl a ) to the eigenbasis {a\ a /3^ CT )- The hamiltonian 
being hcrmitian, the matrix Uka is unitary : UkaU\a — \j\jJha = 1. In the notation U\ a = ^ 



have 



Uk 



\fl + (vo/Vk) 



1 

2 2 



1 + 



we 



(9) 



1 

^2 2 



1 ^fc 
A fc 



The saddle-point equations together with the conservation of the number of conduction electrons are written as 

co = UkVk UF ( E k ) ( 10 ) 



k<7 



1 = ^E^ 



ka 



n c 



kcr 



Their resolution leads to 



|y F | = £>exp[-2/ (fioJc)] 



(11) 



2poo-o/ Ij/f| = 1 



/i = 

where yp = p — ef and po is the bare density of states of conduction electrons (po = 1/2D for a flat band). 
Noting y = E — ep, the density of states at the energy E is p(E) = po (l + (Tq/j/ 2 ). If n c < 1, the chemical 
potential is located just below the upper edge of the a-band. The system is metallic. The density of states at 
the Fermi level is strongly enhanced towards the bare density of states of conduction electrons : p(Ep)/po = 
(1 + cri/yp) ~ l/(2po \vf\)- That corresponds to the flat part of the a-band in Figure 1. It is associated to the 
formation of a Kondo or Abrikosov-Suhl resonance pinned at the Fermi level resulting of the Kondo effect. The 
low-lying excitations are quasiparticles of large effective mass m* as observed in heavy-Fermion systems. Also 
note the presence of a hybridization gap between the a and the (3 bands. The direct gap of value 2(T is much 
larger than the indirect gap equal to 2|?/f|. The saddle-point solution transposes to N=2 the largc-N results 
obtained within the slave-boson mean- field theories (SBMFT). 
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B. Gaussian fluctuations 



We now consider the gaussian fluctuations around the saddle-point solution. Following Read and Newns jl7| , 
we take advantage of the local U(l) gauge transformation of the lagrangian £'(r) 

$i — > Ti exp(«0i) 
fi -> fi exp(i6»j) 
X t -»• A- + i dQi/dr 

We use the radial gauge in which the modulus of both fields $^ and $* are the radial field r.;, and their phase 
(via its time derivative) is incorporated into the Lagrange multiplier A; which turns out to be a field. Use of the 
radial instead of the cartesian gauge bypasses the familiar complications of infrared divergences associated with 
unphysical Goldstone bosons. We let the fields fluctuate away from their saddle-point values : r.; = tq + Sri, 
Xi = Xq + SXi, £f, = <5£ f . and £ Ci = 5£ c . . After integrating out the Grassmann variables in the partition function 
in Equation (0), we get 

Z = J VnVXil)t f .VZ c .exp[-S eff } (12) 

where the effective action is 

k,icj n i i 

with : 

[G -1 (iuj )Y' J ' - ( ^ - * Wn _ ~ ti ^ aa ' _ Jsiifr-Tav'Sij ((T + JcSr^S^'Sij 

[ n)\ ij y ^ o + j c s ri )S aa <5 i:j [-iwn + £f + SXijS^Sij - J s i^ c ..T B 

(2) 

Expanding up to the second order in the Bose fields, one obtains the gaussian corrections S e ^j to the saddle-point 
effective action 

+ ( 8$ 5& ) D^^q, ( |Z \ + ( 6ZJ 5£ ) D^ -1 (q, ( || \ (13) 
where the boson propagators split into the following charge and longitudinal spin parts 

D *(q iu)=( Jc ^~ J c(^2(q>^)+^m(q>^))] ~Jc^M,iuJu)\ , 14 s 
c v V -Jc<Pi(q,iwv) -tp ff (ti,iw v ) J 

r^W-if ■ \ I Js<p\Aq,ie}„) J s [l + Jsvlffaiw*)] 
D'i (q, iuj u ) =| . 11 n . . . „ I, . LJ 



' Js[l + Js<Pf c (<l,iuv)] J|vcc(q, 



and equivalent expression for the transverse spin part D^ _1 (q, iu> u ). The expression of the different bubbles are 
given in the appendix. The charge boson propagator Dc(q, iuj v ) associated to the Kondo effect is equivalent 
to that obtained in the 1/N expansion theories. The originality of the approach is to simultaneously derive 
the spin propagator d| X (q, iu) v ) and Do _1 (q, iu) v ) associated to the spin fluctuation effects. Note that in the 
magnetically-disordered phase, the charge and spin contributions in S e ff are totally decoupled. 
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III. DYNAMICAL SPIN SUSCEPTIBILITY 



Next step is to consider the dynamical spin susceptibility. For that purpose, we study the linear response Mf 
to the source-term — 2S/.B (we consider B colinear to the z-axis). The effect on the partition function expressed 
in Equation (0) is to change the hamiltonian H'j(t) to H'j B (t) 

icrcr' \ a=x,y,z J V ' j i 

(15) 

Introducing the change of variables = £° — iB/Jg, we connect the f magnetization and the ff dynamical spin 
susceptibility to the Hubbard Stratonovich fields £ t. 

IBM , , 

^ = l^ = -(««> + («)(«> ^ 

Using the expression ( [Tjj ) fot the boson propagator D| _1 (q), we get for the longitudinal spin susceptibility 

X^(q, = T2r ii ; , - ii ; J \ w— — 2 || ; , - (17) 



and equivalent expression for the transverse spin susceptibility X//(Qj iu> v ). The diagrammatic representation 
of Equation (JTt]) is reported in Figure 2. The different bubbles y//(q, iu;„), y> cc (q, iw„) and ipf c (q,iuj v ) are 
evaluated from the expressions of the Green's functions 

G//(k, iu n ) = UfeG QQ (k, iu n ) + i^G^k, iw n ) (18) 
G cc (k, iuj n ) = ^G QQ (k, iw„) +ttfcG^(k, iw n ) 

G c /(k,iw n ) = G/ C (k,iw n ) = -u k v k [G aa (k, iu n ) - G (3/3 (k, iuj n )} 

where G QQ (k, iuj n ) and Gpp(\t,iuj n ) are the Green's functions associated to the eigenstates aj^jO >and /3^ CT |0 >. 
In the low frequency limit, one can easily check that the dynamical spin susceptibility may be written as 

, . s Xaa(q,«4/) + X a p(<l,iu v ) 

Xff (q,lU v ) = =2 ; r= ; r (19) 

for both the longitudinal and the transverse parts. 

n F (E-)-n F (E+ +q ) 



ro^y-K-K+q ' -K^k+q) . ,-,+ p - 

Equation ( |T~0| ) constitutes the main result of the paper from which the whole physical discussion on the q- and 
u>- dependence of the dynamical spin susceptibility follows and comparison with experiments is made. 



6 



IV. PHYSICAL DISCUSSION 



From Equation (|19|), one can see that the dynamical spin susceptibility is made of two contributions 

Xintra 

(<l,iu) u ) and Winter (q>^) 

Xff(q.,iu> v ) = Xintra{<\,iu v ) + Xinter (q, iu v ) (20) 

with 

Xintra(q.,lU„) = ^ ; r= ; r (21) 

1 ■ \ X a pi^^u) 
Xinter [<l,lUv) = ^ 7 : \= 7 : V ( 22 ) 

Xintra (q> 2w„)and Xinter (q, respectively represent the renormalized particle- hole pair excitations within the 
lower a band, and from the lower a to the upper [3 band. The latter expression is reminiscent of the behaviour 
proposed by Bernhoeft and Lonzarich |?j to explain the neutron scattering observed in UPt^ with the existence 
of both a "slow" and a "fast" component in x (q, cj)/o; due to spin-orbit coupling. Also in a phenomenological 
way, the same type of feature has been suggested in the duality model developed by Kuramoto and Miyake [^7j . 
To our knowledge, the proposed approach provides the first microscopic derivation from the Kondo lattice model 
of such a behaviour. The bare intraband susceptibility x aa (q ; w ) is weu approximated by a lorentzian 



r (q) 



(23) 



where p aa = x' aa (q> 0) an< 4 To(q) is the relaxation rate of order \yp \ — Tk- This corresponds to the Lindhard 
continuum of the intraband particle-hole pair excitations x aa (<Zi w) 7^ as reported in Figure 3. In the same way, 
we propose to schematize the low- frequency behavior (u << wo(q) of the bare interband susceptibility by 

^ 1 (q,o;)=^(qr 1 (l-^y) (24) 

where p a p — X a (q>0) an d w o(q) is a characteristic frequency-scale of the interband transitions. The value 
of u>o(q) is strongly structure-dependent. In the simple case of a cubic band structure e/. = —2t(cosk x + 
cos fcy+cos k z ) (tight-binding scheme including nearest- neighbor hopping), we find a weakly wavevector dependent 
frequency around q = Q of order of loq = 2 \yp \ / (poJc)- The latter result does not stand for more complicated 
band structures as obtained by de Haas-van Alphen studies pSfl combined with band structure calculations in 
heavy- Fermion compounds. In the following, we will leave u>o(q) as a parameter. Figure 3 reports the continuum 
of interband particle-hole excitations x a p v 7^ 0. Due to the presence of the hybridization gap in the density of 
states, the latter continuum displays a gap equal to 2ao, the value of the direct gap at q = 0, and 2 \uf\, the 
value of the indirect gap at q = Q (close to &f). More precisely, we have shown 

2 

XV' (0, w) - 4p 1° = at 2a < ^ < £ (25) 

lOsJU) 1 — 4<7q 



X a p"(Q,u) = 2po 1 + cj2/(2go)2 at 2\y F \ < u < 2D 

Far from the antiferromagnetic wavevector Q = (jr, n, n), X//(q, to) is dominated by the intraband transitions. 
In the low frequency limit, the frequency dependence of Xmtrai^ w ) can approximate to a lorentzian 

" / \ " /_ \ Xintra (q) -Metros (q) fr,r\ 
X//(q,^) « Xintra (q.w) =W ; " T~TT~ ( 26 ) 

UJ 2 + T intra (q) 

with 
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r intra (q)=r (q)(l-7(q)) 



(27) 



' ( \ - A*" fa) 
XintroW - (l-/(q)) 

^fa) = 4xLa(q,0)x^(q,0). One has: X a Q (0,0) - Paa (0) = p(E F ) and x^(O,0) = Po- The contribu- 
tion expressed in equation ( |2^ ) is consistent with the standard Fermi liquid theory. Note that the product 

r intra 

(q)Xmtra(q) = Paa(q)T (q) is independent of I. 

Oppositely, at the antiferromagnetic wavevector Q, x//(q, w) is driven by the interband contribution and we 
get 

x}/(Q,«) « xL ter (Q,^) = ^ ^„ ter r^ (28) 

^ W max> ) -f- 1 inter 

with 

Wmax = W (l - I) (29) 
Winter = PcWC 1 ~ 7 ) 

where wo, /Qa/3, To and I are the values of Wofa), Pa/3fa) and r (q) and 7(q) at q = Q. The role of the 
interband transitions have already been pointed out in previous works pl| . However while the previous studies 
conclude to the presence of an inelastic peak at finite value of the frequency related to the hybridization gap 
whatever the interaction J is, we emphasize that the renormalization of x Q ^(Q,w) into x««ter(Q,w) leads to 
a noteworthy renormalization of the interband gap. Due to the damping introduced by intraband transitions, 
XinteriQ^) takes a finite value at frequency much smaller than the hybridization gap. Both the relaxation rate 
Winter vanishes and the susceptibility x'inter diverges at the antiferromagnetic transition with again the product 
^interXinter independent of I. Remarkably, the value w max of the maximum of Xj„t er (Q, <j)/u) is at the same 
time pushed to zero. Such a behaviour has been effectively observed in Cei_ x La x Ru 2 Si2 || with a reduction 
of T inter and w max respectively by a factor 4 and 6 when x goes from to 0.075 so when getting closer to the 
magnetic instability occuring at x — 0.08. In order to make the comparison more quantitative, we propose to 
deduce the values of uiq and (1 — I) from the experimental data using the equations d29): uiq = ToTi nte r 



J max 



and (1 — I) = uj^^/ (TqT i n t er ). Table II reports the results starting from the experimental values of ro,ri„ ter , 
w mas (respectively noted Tgs, Tjs, ujq in experimental papers) extracted from the INS results obtained in CeCug 
(ref. ||) and Ce\~ x La x RuiSii at x — and x — 0.075 (ref. ||). The predictions for iv and (1 — I) in these 
compounds seem reasonable. The Stoner enhancement factor (1 — /) decreases in Ce\- x La x RuiSii from x = 
to x — 0.075. (1 — /) of CeCue is intermediate between those two systems. 



V. CONCLUSION 



In this paper, we have set up a new approach of the S — 1/2 Kondo lattice model which enlarges the standard 
1/N expansion theories up on the spin fluctuation effects. The latter effects are proved to be essential for 
the behaviour of the dynamical spin susceptibility near the magnetic phase transition. Our approach provides 
a microscopic derivation of the main features assumed in the phenomenogical models of heavy Fermions as 
the duality model. We predict a two-component behaviour of the dynamical spin susceptibility: a quasielastic 
peak typical of the Fermi liquid excitations, and an inelastic peak at a value w max of the frequency which is 
strongly renormalized due to spin fluctuation effects. Outstandingly well, the frequency of the inelastic peak 
is pushed to zero at the antiferromagnetic transition at the same time as the frequency width vanishes. The 
results have been compared to the Inelastic Neutron Scattering experiment data with reasonable predictions for 
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the Stoner enhancement factor (1 — 1) and the characteristic frequency ujq of the interband contribution to the 
susceptibility. Obviously, more experiments are needed for a systematic test. The issue is important since it 
may have implications for the quantum critical phenomena around the antiferromagnetic critical point. Work 
is currently in progress in that direction and will be presented in a forthcoming paper. We expect the two 
underlined modes to have different effects on the critical behaviour with, on the one hand, the first " intraband" 
mode acting as a paramagnon mode as in the Hertz-Moriya-Millis theory S, and on the other hand, additional 
effects brought by the second "interband" mode. 
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APPENDIX 

The expressions of the different bubbles appearing in the expression of the boson propagators (cf. Eq.|l4|) are 
given here (with i=l, 2, m or ff) 

^(q, iu v ) = ifi(q, iuj v ) + ^(-q, -iw v ) (30) 

kcr,iuj n 



ip 2 ((l,iuJu) = —7j G° CQ (k + q,iuj n + iu;v)G' 7 ffQ (k,iLj n ) 



kcr.iuJri 



cpj f (q,iw u ) = -- ^ G' y ffo (k + ci,iuj n + iu; l/ )G' T ffo (k,iuj n ) 

kcr : iuj n 

<P ] lc(<l, iu ») = -p ^ Co (k+q,w„ + iu v )GZ C0 {\s.,iuj n ) 



■rr.toJ,: 



<p\ c {q.,iUv) = Y G a fca (k + q,iuj n + iuj l/ )G' y fco (k,iuj n ) 



ko.iui 



■a.tuj,. 



ip^.(q,iaj v ) = -~ Y Gl C0 (k + c[,iuj n + iuj u )G l CC0 {k,iuj n ) 

k(T,iuJn 

tpf c (q,iw„) = -- Y G 7 fco (k + q,ico n +iu; u )G l fco (k,iu; n ) 

where G^ Co (k, iw n ), G^ g (k,ioj n ) and GJ CQ (k, iuj n ) are the Green's functions at the saddle-point level obtained 
by inversing the matrix Gg(k, r) defined in Equation (Q). 
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TABLE CAPTIONS 



Table I: Values of the single-site and intersite relaxation rates Tss ar >d Tjs and position of the inelastic peak 
u m ax extracted from the Inelastic Neutron Scattering (INS) measurements performed in CeCue (ref. Q) and 
Ce 1 „ x La x Ru 2 Si 2 at x = and x = 0.075 (ref. §). 

Table II: Predicted values of the characteristic frequency-scale u>o for the interband transitions and the Stoner 
enhancement factor (1—1) from the INS data on Tss, Tis an d ojq (respectively noted To, Ti nter and u> max in 
experimental papers) for the same three compounds as in Table I. Note that Ce\- x La x Ru2Si2 at x = 0.075 
is very close to the antiferromagnetic instability while the Stoner enhancement factor (1 — /) for CeCu^ is 
intermediate between that of the two concentrations x = and x — 0.075 of C e\~ x La x RuiSi2- 
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Table I 





T ss (meV) 


Tjs{meV) 


u max (meV) 


CeCue 


0.42 


0.2 


0.25 


CeRu2Si2 


2.0 


0.75 


1.2 


Ce\- x La x Ru2Si2 


1.4 


0.2 


0.2 



Table II 





To(meV) 


r int er(meV) 


u m ax{meV) 


wolded.Cmey) 


(l-/)|«fed. 




0.42 


0.2 


0.25 


0.34 


0.74 


CeRu 2 Si2 


2.0 


0.75 


1.2 


1.25 


0.96 


Cei- x La x Ru2Si2 


1.4 


0.2 


0.2 


1.4 


0.14 
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FIGURE CAPTIONS 



Figure 1: Energy versus wave- vector k for the two bands a and /?. Note the presence of a direct gap of value 
2<7o and of an indirect gap of value 2 \uf\- 

Figure 2: Diagrammatic representation of Equation (poj ) for the dynamical spin susceptibility Xffili uj )- 

Figure 3: Continuum of the intra- and interband electron-hole pair excitations \ aa {q> w ) 7^ and x (q, lo) ^ 0. 
Note the presence of a gap in the interband transitions equal to the indirect gap of value 2|yi?| at q = kp, and 
to the direct gap of value 2ao at q = 0. 
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FIG. 2. 
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FIG. 3. 
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